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Abstract 

The results of our study of the motion of a three particle, self-gravitating system in general relativistic 
lineal gravity is presented for an arbitrary ratio of the particle masses. We derive a canonical expression 
for the Hamiltonian of the system and discuss the numerical solution of the resulting equations of mo- 
tion. This solution is compared to the corresponding non-relativistic and post-Newtonian approximation 
solutions so that the dynamics of the fully relativistic system can be interpretted as a correction to the 
one-dimensional Newtonian self-gravitating system. We find that the structure of the phase space of each 
of these systems yields a large variety of interesting dynamics that can be divided into three distinct 
regions: annulus, pretzel, and chaotic; the first two being regions of quasi-periodicity while the latter is 
a region of chaos. By changing the relative masses of the three particles we find that the relative sizes 
of these three phase space regions changes and that this deformation can be interpreted physically in 
terms of the gravitational interactions of the particles. Furthermore, we find that many of the interesting 
characteristics found in the case where all of the particles share the same mass also appears in our more 
general study. We find that there are additional regions of chaos in the unequal mass system which are 
not present in the equal mass case. We compare these results to those found in similar systems. 
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1 Introduction 



The calculation of the motion of N particles under a specified, mutual force is one of the oldest problems 
in physics, commonly referred to as the iV-body problem. This problem occurs frequently in many distinct 
subfields and remains an active area of research. When the specified force is that of Newtonian gravity in 
three spatial dimensions, a closed form solution of the motion can be obtained for N = 2. This is not true 
for (3+1) dimensional general relativistic gravity, however, due to the existence of energy dissipation in the 
form of gravitational radiation. All attempts to calculate the motion of more than one particle in (3+1) 
general relativity have required some form of approximation. 

Considerable progress in this area of research has been made recently by reducing the number of spatial 
dimensions from three to one. These lower-dimensional theories provide a simpler prototype for their higher- 
dimensional counterparts. Furthermore, for Newtonian gravity, one-dimensional self-gravitating systems 
(OGS's) have proven to be very useful in modelling many diverse physical systems. For example, it has been 
found that there exist stable core-halo structures in the OGS phase-space that are reminiscent of those found 
in globular clusters T . These structures consist of a dense core of particles near equilibrium surrounded by 
a cloud of high kinetic energy particles that interact very weakly with the core. The OGS also models the 
dynamics of flat, parallel sheets colliding along a perpendicular axis |2] and the motion of stars interacting 
with a highly flattened galaxy 0. More specifically, the three-particle OGS has been found to model the 
motion of a billiard colliding with a wedge in a uniform, gravitational field j2 , two elastically colliding billiard 
balls in a uniform, gravitational field U, and a bound state of three quarks to form a "linear baryon" [5]. 
There are still many open questions about the OGS concerning its ergodic behaviour, the conditions (if any) 
under which equipartition of energy is attained, whether or not it can reach a true equilibrium configuration 
from arbitrary initial conditions, and the appearance of fractal behaviour j2J. 

In a relativistic context, reduction of the number of spatial dimensions results in an absence of gravita- 
tional radiation while retaining most (if not all) of the remaining conceptual features of relativistic gravity. 
Consequently, one might hope to obtain insight into the nature of relativistic dynamical gravitational sys- 
tems at the classical (and perhaps even quantum) level in a wide variety of physical situations by studying 
the relativistic OGS, or ROGS. 

Comparatively little has been known about the ROGS (even for N = 2) until quite recently, when a pre- 
scription for obtaining its Hamiltonian from a generally covariant, minimally-coupled action was obtained 
j7|. In the non-relativistic limit (c — *■ oo), the Hamiltonian reduces to that of the OGS. This opened up 
the possibility of extending the insights of the OGS into the relativistic regime, and indeed, considerable 
progress has been made. Exact, closed-form solutions to the 2-body problem have been obtained jH]- These 
have been extended to include both a cosmological constant EI| and electromagnetic interactions jll| . 
and a new exact solution to the static-balance problem has been obtained In the iV-body case, the 

Hamiltonian can be obtained as a series expansion in inverse powers of the speed of light c to arbitrary order 
and a complete derivation of the partition and single-particle distribution functions has been found in both 
the canonical and microcanonical ensembles |13| providing interesting information concerning the influence 
of relativistic effects on self-gravitating systems. Very recently, formulation of the ROGS has been extended 
to circular topologies ^A^ (forbidden for the OGS), and a new A^-body dynamical equilibrium solution has 
been found jJS]. An exact expression for the relativistic, 3-body Hamiltonian has been calculated and the 
motion of three equal-mass particles has been extensively studied [161 117| : these results are summarized in 

§o 

In this paper, we will generalize the study of the motion of three particles to include the unequal mass 
case. We work with a 2D theory of gravity on a line (lineal gravity) that models 4D general relativity in 
that it sets the Ricci scalar R equal to the trace of the stress-energy of prescribed matter fields and sources. 
Hence, as in (3 + 1) dimensions, the evolution of spacetime curvature is governed by the matter distribution, 
which in turn is governed by the dynamics of spacetime ^H]- Sometimes referred to as i? = T theory, it is a 
particular member of a broad class of dilaton gravity theories formulated on a line. What singles it out for 
consideration is its consistent non-rclativistic (i.e. c — > oo ) limit |18| which is, in general, a problematic limit 
for a generic (1 + l)-dimensional theory of gravity jl9| . Consequently it contains each of the aforementioned 
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non-relativistic self-gravitating systems as special cases. Furthermore, it reduces to Jackiw-Teitelboim (JT) 
theory [201 when the stress-energy is that of a cosmological constant. 

We have found that the best way to study the motion of three particles is to work in the canonical 
formalism. By expressing the action in canonical variables we are able to determine the Hamiltonian as 
a spatial integral of the second derivative of the dilaton field. This field is determined by the constraint 
equations derived from the action which can be solved by matching the solution of the field across each of 
the three particles. The result is a transcendental equation containing the Hamiltonian and expressed in 
terms of the remaining degrees of freedom, that is, the two mutual separations of the particles and their 
conjugate momenta. From this transcendental equation we obtain the canonical equations of motion, which 
are then solved numerically. 

Through a change of coordinates, the Newtonian, three-particle OGS can be shown to be isomorphic to 
the motion of a single particle in a linear, hexagonal well potential. By applying this same change of variables 
to the three-particle ROGS we find an analogous hexagonal potential where the sides of the hexagonal cross- 
section are curved outwards and the sides of the well no longer increase linearly with increasing particle 
separation. We find that, by changing the relative masses of the particles, the shape of the hexagonal cross- 
section in both the Newtonian and relativistic systems is expanded or contracted perpendicular to one of the 
lines connecting opposite vertices. This change of variables simplifies the analysis of the motion significantly 
and is used throughout to extract useful information from both the three-particle OGS and ROGS. 

As in ref. ^7] we consider three distinct 3-body, self-gravitating systems: the non-relativistic case (N) 
which has been extensively studied in many different contexts jHEEH^^, the fully relativistic case (R) 
described above, and a post-Newtonian expansion (pN) of the R system, truncated to leading order in c~^. 
While exact relativistic solutions of the N-body problem have only been found for N = 2,3, the post- 
Newtonian expansion has been found for all finite values of N up to any order of accuracy [7] . Both the R 
and pN systems reduce to the N system in the limit c — *■ cx). 

In section|21we outline the canonical reduction procedure of ^21 that leads to the relativistic Hamiltonian 
expression and the resulting canonical equations of motion. Some general properties of each of the systems 
are then discussed in section |3 focussing on the character of the associated potential energy functions of 
each. The method for numerically solving the equations of motion is described in section^jand the results of 
this numerical solution presented in section|Sl These results are then summarized and discussed in sectional 
which concludes with a comment on areas of further research interest. 



2 Hamiltonian Formulation of the Relativistic Equations of Mo- 
tion 

The general procedure for deriving the iV-body Hamiltonian via canonical reduction is given in [HI 1101 117j 
and only a brief description will be given here. 

The action for the gravitational field minimally coupled to N point particles in (1+1) dimensions is given 

by 



N 
a=l 



(1) 



dz^ dz 

a, I 
dTa dTa 




where g^^i, is the metric tensor with determinant g, R^i, is the Ricci tensor, Tq the proper time for the 
ath particle with mass uia and position Za, and k = SttG/c'*. We use to denote the covariant derivative 
associated with g^^,^. The scalar (dilaton) field has been incorporated because the classical Einstein-Hilbert 

^ These studies examine the 3-body problem in a classical potential obtained by solving Poisson's equation in one spatial 
dimension. This potential linearly depends on the separation of the particles as seen in 1201 . The chaotic properties of the 
one dimensional 3-body problem with a potential that depends inversely on the separation (as in three dimensions) have been 
studied in 
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action in (1+1) dimensions is trivial due to the vanishing of the Einstein tensor. This action describes a 
self-gravitating system of N particles without coUisional terms {i.e. the particles pass through each other). 
^From Q one can derive the following field equations 

i?-5^"V,,V,* = 0, (2) 



d ( , .dzn 1 , ^dz'^dz^^ 



a ^ a 



dTa I dTa ) 2 dTa dTa 



0, (4) 



where 



Tf_,u ^y^ma fdTa—^g^,agwp'^T^'^5^^Hx~Za{Ta)) (5) 
„ J dTa dTa 



is the stress-energy tensor for the N particles and is conserved via Inserting the trace of l|31l into Q we 
obtain 

R^kT'^^. (6) 

The fact that we retain this simple relation between the geometry of spacetime and the matter, analogous 
to the Einstein field equations, is the motivation for choosing the dilaton coupling in (|T|). 

Equations l@J and © form a iV -f 1 system that can be solved for the single metric degree of freedom and 
the TV particle degrees of freedom. Equation Q relates the evolution of the dilaton field to the evolution of 
the point masses. 

To arrive at a Hamiltonian theory we begin by writing the metric as 

ds^ = -N^ (x, t) df +-i{dx+ ^dtj (7) 

where Nq and A^i are the lapse and shift functions which act as Lagrange multipliers for the resulting 
constraints of the Hamiltonian system and 7 is the single metric degree of freedom. 

By also defining pa, tt and 11 to be the conjugate momentum of Za, 7 and 4" respectively, one can 
canonically reduce the action to the form [Jj 



I = J d^x i^y^PaZaSj 



x-z„) + iA*l (8) 



K 



upon eliminating the constraints and choosing the coordinate conditions 7=1 and 11 = 0. Here we use A 
to denote d^/dx^ and a dot to denote d/dt. With the action in this form, we recognize the second term 
7i = — ■i-A^E' as the Hamiltonian density and can immediately write down the Hamiltonian for N particles 
as 

[ dxH = -- / dxA^ (9) 



where is a function of Za and Pa and can be determined from the solution to the constraint equations 
which now take the form 

A* - ^{^'f + K^Tl^ +>i^\/pl + ml5{x - Za) = (10) 

a 

+ Y,Va5{x-Za)^Q (11) 

a 

where a prime denotes d/dx. 

The solution of , Hl()|l , and for the 3 particle case is given in Tf and will not be reproduced here in 
detail. The basic procedure involves choosing a specific configuration of the 3 particles and solving ((Tn|l and 
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pifl in the region between each particle. The constants of integration are then determined by demanding 
that ^ and 4"' remain finite and coincide at the position of the particles. This gives an implicit equation for 
the Hamiltonian H for the specified particle configuration and the Hamiltonian for a general configuration 
is obtained by permutation of the particle indices (1, 2, and 3). 
This implicit equation for the Hamiltonian can be expressed as 

L1L2LS = Mi2M2iL*3e^'^^^^^^+-^^'^^^''-^^^+^^^'>^^^^ 

+ 7\/t23Al32L*e^^^3[(i2+Al23)z21-(i3 + A^32)23l] 

+ 7W3i7Wi3L2e*""''^^'+'^'''^''"^^'+-^"^^''l (12) 

or more compactly 

ijk 

where 

M.,j ^ Mr - epis.j, M, = \jvi + ml (14) 

U = H-M,-ei^V,s,,) L* = {\- W s,,s,k)M, + Li (15) 

j j<k=ii 

with Zij = {zi Zj) , Sij — sgn{zij), and e'-''^ is the 3-dimensional Levi-Civita tensor. The discrete parameter 
e = ±1 is a constant of integration that flips sign under time reversal. This provides a measure of the flow 
of time of the gravitational field relative to the particle momenta. 

Although we cannot obtain an explicit expression for the Hamiltonian, we are able to derive the equations 
of motion explicitly by partially differentiating (|13|l implicitly with respect to Za and Pa and solving for 
dH/dza and dH/dpa respectively. From Hamilton's equations 

dH , , 

ia = ^ 16 

OPa 

OZa 

we can obtain the equations of motion. 

For example, for a = 1, and p7|l become 

ii I L2LZ + LiLz + L1L2 

- [M2 - ep2S2i][M, - episi2][l + jL*|^i2|]et--[(^i+-^-)-i3-(L2+A^2i).23] 

- [M3 - ep3S3i][Afi - episi3][l + ^L*|^i3|]eT-i3[(Li+M.3).i2+(L3+>(32).23] 

-[M2 - ep2S2s]Wh - ep3S32][l + jL*|z23|]et^=3[(L3+X3i).i3-(L2+Mi3).i2]| 

= [Ah - ep2S2i][(^ - esu)Ll - (Ah ~ episi2)(esi3 + |L*(ezi2))]e^--[(ii+-W-)-i3-(L.+.M.i)^23] 
+ [M3 - eP3S3i][(^ - esi3)L* - {Ah - epi5i3){esi2 + jL*(ezi3)}]e^^i3[(Li+.M.3).i2+(L3+.M32).23] 

+ [M2 - ep2S23][M3 - ep3S32][-Sl2Si3^ + \s23L\\^\zr2\ ~ £ |zi3 1]] ef 1(^3+^31 ).13 -(L.+>( 13)^12] 

+ ^^1.2^3 + e(si2LlL3 + 513^2-^^1) (18) 
dpi 
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and 



T 

Pi 1 L2 


L3 + LIL3 + 


LIL2 






- [M2 


1 r 71 r 

- eP2S2l\[Mi 


- episi2j[l 


+ 5^3^121 


Jgf S12[(il+A1l2)zi3-(i2+A^2l)z23] 


- m 


- ep3S3i][Mi 


- episi3][l 


+ ^L*2\zi3\ 


]gf Sl3[(il+A^23)zi2 + (i3+-M32)z23] 


-[M2- 




- eP3S32][l - 


^^LI\Z2S\ 


J gK/4s23[(i3 + -fl-f3-eP3S32) 213 -(i2+Af2-ep2 823)212] 1, 


= [M2 


- ep2S2i][Mi 


- episi2j[^ 


si2Ll[H + 


- e(p2 - Pi)si2 + ep3Si3]]eT''i2[(^i+-^i2)^i3-(L2+Al2i)223 




- ep3S3i][Mi 


-episi3][^ 


sizLl[H + 




+ [M2 


- ep2S23][A^3 


- eP3S32][^ 


S2iL*lPl{si2 - Si3)]eT«23[(i3+A^3i)^i3-(L2+A1i3)^i2]^ 



(19) 

The equations for a = 2, 3 are similar and will be omitted here. 

3 General Properties of the 3-Body System 

Before we go on to solve the equations of motion, it is instructive to consider some general characteristics 
of the 3-body system described by the determining equation p2|l and its associated non-relativistic and 
post-Newtonian counterparts. 

To compare the relativistic motion to that predicted classically, we introduce the Newtonian N particle 
Hamiltonian in (1+1) dimensions 

a a b 

where Zab = Za — Zb as before. We shall refer to this as the Newtonian or N system. 

A post-Newtonian approximation of the general A^-body Hamiltonian has been found "7" and is given 
here up to order 

HpN = X! "^"^ + X! 2^ ^ ~ ^ ^ "^a^h \Zab\ H ^ ^ ^ {maPb ~ mbPa) {Zab) ^ 

a a ah ah 



a b 

4 V4 



^ ^ ^ mambnic {\Zab\ \Zac\ + ZabZac) ■ (21) 
a b c 

If we re-scale ()21f) to remove the constant rUaC^ term and take the limit c — > cx) then it is clear that we 
retrieve the Newtonian result H2U|) . 

However, the coordinates Za and Pa are not necessarily the most natural coordinates to use to describe 
the post-Newtonian system. The reason is that the fourth term on the right hand side of H21I) is proportional 
to c~^. In (3 + 1) dimensions, terms in odd powers of c^^ are associated with gravitational radiation, but, 
in (1 -I- 1) dimensions, there are not enough degrees of freedom to allow for the existence of gravitational 
radiation, and such terms are artifacts of the choice of canonical variables. 

Indeed, as in jTj, we can remove the c^^ term by performing the canonical transformation 

Za ^ Za Za (^^) 



en \ ^ 

Pa ^ Pa ^ Pa - ^ ^ maTTlbZab 
b 



(23) 
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after which, (|21|l becomes 



H, 



pN 



Pi 



Pi 



KC 



EE 



8to3 



^ 2rna 

a 



JTlaTOfc Fob I 



nib ■ 



E E E "^ambrric i\zab\ 



(24) 



where = — z^. Since H24|) uses different coordinates than (|20|l and p2|) it is important to distinguish 
between the two expressions for the post-Newtonian Hamiltonian. We wiU refer to the system described by 
H21|) as the untransformed post-Newtonian system, or UpN system, and the system described by H24|l simply 
as the post-Newtonian or pN system. Note that only the pN system was studied in ^7]. To complete our 
nomenclature, the fully relativistic system will be denoted as the R system. Unless otherwise stated, the 
Newtonian system will be assumed to have been re-scaled so that H{zab — 0,pa ~ 0) ~ {mi + m2 + 7713)0^ 
in order to properly compare it to the pN and R cases. 

In order to simplify our analysis we will adopt the convention of 5 and 17 and define the following 
canonical coordinates: 



1 



p = —1= {zi - Z2) , X 



V2 



1 , 



Zl - 2^3) , Z = Zi -I- Z2 -I- Z3, 



with conjugate momenta 
1 



1 



Vi) 



PA = [p\ +P2- 2p3 



= g {Pi 



P2 +P3) 



(25) 



(26) 



In the non-relativistic limit, Z and pz are related to the center of mass and its conjugate momentum. While 
the equivalence principle does not allow us to arbitrarily set Z in the relativistic case, we can, without loss 
of generality, choose pz to vanish. The consequence of this is that we can explicitly express pi, p2, and p^ in 
terms of the newly defined momenta (|26(l but can only express the separations of the particles Zab explicitly 
in terms of the new coordinates H25(l . This gives us the following relations: 



2^12 



Zl3 



11 11 2 

Pp, P2 = —^P\ l=Pp, P3 = -\l 7;PX- 



(27) 



(28) 



72""' V6" %/2""' V3^ 

All of the Hamiltonian expressions (|12|1 . (|2U|) . H21I) . and H24|) do not depend on Z or pz and so these 
variables are irrelevant. Expressions for (|12l) . (|20|l . and H24|l in terms of the new coordinates are given in [T7| 
for the case when mi = m2 = m^. The corresponding expressions for unequal masses are very cumbersome 
and will not be reproduced here. 

By defining the potential of each system as V{p,X) = H[pp = 0,pa = 0) we can compare some of the 
different characteristics of the 3 systems. In p-X coordinates the N potential becomes 

Vn = mi +m2 + ms H — ( 2mim2 \p\ + mim^ 



4V2 



Vsx 



m2mz 



V3A 



P 



(29) 



where we have rescaled the Hamiltonian as described above (with c henceforth set to unity unless explicitly 
stated otherwise). The UpN potential is given as 



VpN = mi + m2 + m3 + 



4V2 V 



( 2r7iim2 IpI -l- TOim3 -y/S A + 



P 



+ 77127713 



V3A 



^ (l) ^ ™iTO2TO3 (^V + (n/3A % (V3A - V (l-f Spsi) \p\ VSA -I- 
^2(4) ™i"^2™3 ( (1 - SpS2) IpI VSA - p + - (1 + S1S2) VSA 4- /9 VSA- 



(30) 
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and the pN potential is 



VpN = mi + 7712 + ma + 



4\/2 



2TOim2 IpI + mima 



+ 17121713 



1 /K\2 



2 V4 



717117121713 (^{1 - SpSl) \p\ 



+ {l + SpS2)\p\ 



V3A 



(31) 



where p and A are defined as in (|25|l using the Za coordinates of H22|l . Here we have defined Sp = sgn{p), 
si — sgn (V3A + p) , and S2 — sgn (-y/SA — p) and the s terms are defined similarly in terms of p and A. The 
exact rclativistic potential can be calculated from p2|l to be 



{Vr - mi) {Vr - TO2) {Vr - TO3) = TO1TO2 (Vr - siS27n3) exp 

+ mima {Vr + SpS27n2) exp 
+ m2m3 {Vr - SpSimi) exp 



.2^/2 

K 

4^2 

K 



Vr\p 

Vr 



Vr 



Vsx + p 
Vsx- p 



(32) 



An extensive comparison between the different potentials has been given in |17| for the case where the 
particle masses are equal and so here we wish to focus on the changes to the potential due to changes in the 
relative masses of the 3 particles. 

A cross section of each of the potentials at a fixed value of V is shown in Figure for the case where all 
particles have the same mass. All of the potentials share a certain hexagonal symmetry in that they are all 
smooth except along the lines 



p = 
p + V3A = 
p - Vsx = 



(33) 
(34) 
(35) 



which correspond to zi — Z2, zi — Z3, and Z2 = Z3 respectively {i.e. the potential is not differentiable when 
two particles are coincident). This is true for all ratios of the masses of the particles. 

The Newtonian potential is a distorted hexagonal well with sides that increase linearly with V. The 
hexagonal cross-section at any value of Vn only has equal length sides when mi = TO2 = ms. Figure |21 shows 
various cross-sections of the Newtonian potential at a fixed value of V for different mass ratios. We see that 
increasing the mass of particle 3 has the effect of expanding the the hexagon away from the p = bisector 
while decreasing the mass contracts the hexagonal cross-section towards p — 0. Increasing and decreasing 
the mass of particles 1 or 2 has the same effect but the deformation is perpendicular to the p — \/3X = or 
p + a/SA = respectively. When all three particles have unequal mass, the hexagon is deformed as above 
with the magnitude of the deformation in each direction given by the relative values of the mass. 

The rclativistic potential is similar to the Newtonian potential except that the sides of the hexagon 
become concave. Furthermore, for small values of (p, A), the rclativistic potential increases much more 
rapidly than the Newtonian potential. However, at a value Vr such that 




- mjj I^Vr - (A/tot - rrij) ^ 
(Aftot - TTij) rrij 



Vr 



1 



{Vr 



Vr - (Aftot - mj) 



(36) 



(for j = 1,2 or 3), the slope of the rclativistic potential becomes infinite, after which the size of the 
distorted hexagon decreases like (IuVr) /Vr with increasing Vr. In the equal mass case this yields a value 
Vr ~ 6.71197mc^, where m = Aftot/3. For m = Aftot/2 we obtain k, 6.886682mc^ which is the maximal 
possible critical value of the potential, and in the limits m — > 0, Aftot we find Vr — > Aftot- A plot of the 
critical values of the potential as a function of m^is given in fig. O 
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p 



Figure 1: A cross section of the 4 potentials at V « 1.3MtotC^ in the case that all 3 particles have the same 
mass. N - solid, pN - dashed, UpN - dotted, R - dash-dotted. The p and A are dimensionless variables 
defined using the dimensionless positions £j of equation (|38|l . 




-1.5- ^ 



-2 - 

-2.5 1 ' ' J 1 ' 

-3-2-10123 

P 

Figure 2: Cross-sections of the Newtonian potential at « l.SMjotc^ for various mass ratios 1711:1712:1713. 
Solid - 1:1:1; dashed - 1:1:4; dotted - 4:4:1; dash-dotted - 1:4:8. Notice that all discontinuities lie on one of 
the three bisectors or regardless of the mass ratio. The p and A are dimensionless variables 

as in Figure n 
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V 
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0.5 1 1.5 2 2.5 3 

m 



Figure 3: Critical values of the relativistic potential Vr as a function of a given particle mass in units of 
Mtot (here set equal to 3). The maximum critical value occurs in the case when nij = Mtot/^- The minimal 
value approaches the limit Vr — M^ot as mj — > or Mtot- 
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-1 -0.5 0.5 1 

P 

Figure 4: Cross-sections of the relativistic potential at l.^MtotC^ for different ratios of the particle 

masses. The correspondence between line style and ratio is the same as in Figure|21 The deformation of the 
potential due to changing the mass ratio is the same as in the Newtonian case. The p and A are dimensionless 
variables as in Figure Q 

The overall shape of the relativistic potential is that of a distorted, hexagonal carafe. The distortion of 
the relativistic potential for different ratios of the particle masses is analagous to the Newtonian potential 
and can be seen in Figure 01 

The untransformed post-Newtonian potential shares similar features with the relativistic potential in 
that the sides of the distorted hexagonal cross-section are concave outward. However, as one might expect, 
the potential increases less rapidly than the relativistic potential but still more rapidly than the Newtonian 
potential at small values of (p, A). Furthermore, the sides of the well continue to increase quadratically 
with increasing (p, A) without the slope ever going to infinity as in the relativistic case. Figure |5l shows a 
cross-section of the untransformed post-Newtonian potential at a fixed value of V for different mass ratios. 

The transformed post-Newtonian potential has a much different character than all of the potential energy 
functions discussed so far. The sides of the distorted hexagon become convex and the vertices are always 
coincident with the Newtonian potential at a fixed value of V. As V increases, the sides become more convex 
with respect to the Newtonian potential. Cross-sections of this potential for different ratios of the particle 
masses at a fixed value of V can be seen in Figure El 

Finally, we note that the potential energy does not completely govern the motion in the R and pN cases 
as it does in the N case due to the momentum dependence of V in the former cases. Consequently, such 
comparison of the potentials is limited in the insight it can provide. 

4 Methods for Solving the Equations of Motion 

The motion of the 3 particles under study is quite complex so we have adopted several methods to study 
the equations of motion. The most straight forward approach is to look at the position of each particle with 
respect to the center of mass J2 

function of time (where time will be explicitly defined shortly) . 
Recall, however, that the choice of the center of mass reference frame is arbitrary and is not necessarily 
stationary to an observer as it is in the Newtonian case. 
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Figure 5: Cross-sections of both the untransformed (left) and transformed post-Newtonian potentials at 
V « l.ZMtot^? for various mass ratios. The correspondence between ratio and line is the same as in 
Figure [3 The p and A are dimensionless variables as in Figure ^ 

Under the change of coordinates H25|l and (|26() , the motion of the 3 particles is isomorphic to the motion 
of a single particle moving in the hexagonal well potential (|29|l in the Newtonian case. The situation is 
analogous in the fully relativistic and post-Newtonian cases except that the potentials become momentum 
dependent. So, as in ^7] and [S], we focus on the trajectory of this particle, which we call the 'hex-particle', 
in the p-A plane as an alternate way to analyze the motion. 

As mentioned before, the bisectors correspond to points where two of the particles are coinci- 

dent. These bisectors divide the p-A plane into 6 sextants corresponding to the 6 different configurations the 
3 particles can assume. So, when the hex-particle moves from one sextant to the next, this corresponds to two 
particles passing through each other. In the equal mass case, all 6 sextants are equivalent, and in the unequal 
mass case, opposite sextants correspond to the opposite configuration of particles {i.e. (1, 2, 3) (3, 2, 1)). 
Further symmetries exist when two particles have the same mass. 

An analogous system was studied by Lehtihet and Miller |2 1^ who demonstrated that a self-gravitating 
Newtonian system of 3 particles in (1-1-1) dimensions with collisions is equivalent to the motion of a particle 
in a uniform gravitational field colliding elastically with a wedge. The existence of particle collisions in 
our study of the N system would correspond to the hex-particle being confined to a single sextant where it 
would reduce to the particle-wedge system. In this particle-wedge system, the equations of motion can be 
integrated between collisions of the particle with the wedge and a discrete mapping describing the radial 
and angular velocity of the particle at each collision can be used to describe the motion. The simplification 
to a discrete mapping allows one to calculate fixed points in phase space and evaluate their stability much 
more easily. Unfortunately, the equations of motion for the pN and R systems are much more complex than 
in the Newtonian case and it is not clear how to create a discrete mapping between particle collisions. 

Following and |17| we define two types of motion: A motion, where two particles cross twice in 
succession; and B motion, where a single particle crosses each of the other two in succession. In (p. A) 
space, after the hex-particle has just crossed one of the bisectors, A motion would correspond to a crossing 
of the same bisector while B motion would correspond to it crossing a different bisector. In this way, one 
can describe the trajectory of the hex-particle as a succession of A and B motions and develop a 'symbol 
sequence' for a given trajectory. To simplify the notation, we use exponents to denote a number of repeats of 
a given type of motion so that the symbol sequence takes the form fl- ■ ^ )''' where hj^rii, rij S 

and Ik is possibly infinite, in which case we denote it by an overbar {i.e. linic^oo{A'^ B^^ = A°-B''). Since 
the type of hex-particle motion at a given bisector depends on the previous bisector, we avoid ambiguity by 
saying that the first bisector crossing of the hex-particle is undefined, and the symbol sequence begins at 
the second crossing. To aid in understanding this nomenclature we have listed the symbol sequence in the 
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captions of all configuration space trajectories where the trajectory is easy to follow. 

The above methods allow us to study and classify individual trajectories of the hex-particle in the 4 
dimensional (p, X,pp,px) phase space. In order to study some of the global structure of this phase space, we 
construct Poincare maps. Since all of the Hamiltonians under study are time independent, the total energy 
of the system is a constant of motion and so the motion at a given energy is confined to a 3 dimensional 
hypersurface in phase space. We can further reduce this to 2 dimensions by plotting the radial momentum 
PR and the square of the angular momentum pg of the hex-particle each time it crosses one of the bisectors, 
as in [2] and ^7]. This is known as the surface of section, or Poincare map. 

In the equal mass case, all bisectors are equivalent and so pn and pg at each bisector can be plotted on the 
same surface of section, as in . When the masses are unequal, this procedure is not possible and one must 
distinguish between the different bisectors and the directions in which they cross. We have chosen to plot 
points on the Poincare map each time the hex-particle crosses the p = boundary in the positive angular 
direction (i.e. when pg > 0). In the case when mi — m2, particles 1 and 2 are indistinguishable and we 
may also plot {pii,pg) each time the hex-particle crosses the p = bisector in the negative angular direction 
{pe < 0). Due to the nature of the Hamiltonian phase space, the different surfaces of section corresponding 
to the different bisectors and directions contain the same information and no generality is lost in making the 
above choice. 

Since we were unable to find a closed form solution to either the relativistic determining equation H12() or 
the equations of motion and (|17|1 . it was necessary to employ numerical techniques to study the motion. 
Using a Matlab integration routine (ode 15s) we were able to solve the equations of motion in the N, pN, 
and R systems. 

The odel5s routine uses a variable order method for solving stiff differential equations 23 . In order to 
control computational errors, we imposed absolute and relative error tolerances in the numerical routine of 
Cahs = irei = 10^^ SO that the estimated error in each of the dynamical variables p{i), X{i),pp{i), and p\{i) 
at each step i in the numerical integration is 

e{i) <uia-x{erei\yii)\,eabs) (37) 

where y{i) represents a generic, dimensionless dynamical variable at time step i. These dimensionless vari- 
ables will be introduced shortly. 

Furthermore, we periodically checked that the total energy of the system remained constant to ensure 
that the solution was stable and physically correct. 

We found (both in this study and in ref. ^7]) that the numerical precision available to the computer did 
not allow the integration routine to solve the equations at energies approximately H > 2Mtot(? ■ We were 
unable to find a numerical integration routine that could integrate the equations of motion in this energy 
regime so the dynamics of the system at high energies still remains an open problem. 

Furthermore, when we integrate the pN equations of motion, we find that the resulting energy of the 
numerical solution does not remain constant in time, despite the fact that the Hamiltonian H24|l describes 
a conservative system. The variation in energy becomes greater as the differences between the masses 
increases. For the case when all masses are equal, this variation is on the order of the imposed numerical 
error tolerances and can be ignored. A description of the dynamics of the pN system in the equal mass case 
is given in |17j . The variation in energy increases drastically when we change the ratio of the masses even by 
a small amount. For example, when we integrate the equations of motion in the case where the mass of one 
particle is half that of the other two, we see a variation in the total energy on the order of 10~'^MtotC^ over 
the duration of the trajectory. The cause of this energy variation is unclear but its magnitude is clearly too 
large to ignore. Due to this energy fluctuation, the numerical solutions to the post-Newtonian equations of 
motion that we obtained in the unequal mass case are clearly unphysical and will not be presented in this 
paper. 

We cast the expressions for the Hamiltonian and the equations of motion in the different systems in a 
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dimensionless form using the coordinates Zi and pi, given by 

4 



(38) 

Pi = Mtotcpi- (39) 
We then express the dimensionless Hamihonian as 

Mtotc^ 

so that 77 corresponds to H being equal to the total rest mass of the system. The total, dimensionless 
energy for all systems is ry + 1 (recall that we are assuming the Newtonian Hamiltonian has been rescaled 
so that the zero point is the total rest-energy of the system) . In this way, a single value of rj corresponds to 
the same energy in all 3 systems. 

The equations of motion then become 



dr] IdH _ 4 dzi _ dzi 

dpi c dpi nMtotC^ dt di 

drj 4 dH 4 dpi dpi 

dzi kM^^^c'*' dzi nMtotc^ dt di 

where we recognize t as the dimensionless time unit, given as 



(41) 
(42) 



t = J - (43) 

We refer to < = 1 as one time step, p. A, Pp, and p\, the dimensionless counterparts of p, A, Pp, and px 
respectively, are defined as in (^5]! and using the hatted variables of lp!H|l and In the subsequent 

analysis, dimensionless variables will be assumed unless otherwise stated. 



5 Solution to the Equations of Motion 

In this section we present the results of our numerical analysis of the equations of motion. In § 15.11 we 
summarize the equal mass results of | 17| then go on to present how the dynamics change in the unequal 
mass case in § 15.21 and § 15.31 

5.1 Equal Mass Solution 

The study of the N, pN, and R systems when all three particles share the same mass revealed a large variety 
of different types of trajectories. The different types of motion can be classified into 3 broad categories 
which we call annulus, pretzel, and chaotic. Note that our naming scheme is not standard in the literature 
of dynamical systems. Our nomenclature was chosen because of its direct physical interpretation in terms 
of the three particles (despite the fact that the terms annulus and pretzel derive from the shape of the 
trajectories in the rho-lambda plane). 

Annulus trajectories correspond to the hex-particle never crossing the same bisector twice in a row 
resulting in an orbit about the origin of the {p, A) plane. The symbol sequence for all annulus orbits is B. 
In terms of particles, these trajectories represent all motions in which no two particles cross successively. 
Most of the trajectories in this class never exactly repeat themselves after any number of orbits about the 
origin. The result is a densely filled region of the (p. A) plane circling the origin. All of these trajectories 
form closed loops on a Poincare map. 

Pretzel trajectories are so named because of the complex patterns they make when plotted in (p. A) 
coordinates. Symbolically, these trajectories fall into two distinct classes: 1) regular trajectories, which are 
denoted by some repeating pattern of A's and S's {e.g. B^"^ A'^ B^) and 2) quasi-regular trajectories, 
represented by some repeating pattern of A's and B's with extra A motions occasionally occurring on each 
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repetition of the pattern {e.g. A'^ {A^ B^)^ {A^ B^y^ . . .). As in the annulus case, most pretzel trajectories 
never exactly repeat themselves and densely fill a region of phase space. Pretzel trajectories appear either 
as a series of small enclosed loops or as a series of disconnected lines on a Poincare map. 

Chaotic trajectories are those that eventually cover all allowed regions of phase space and are denoted 
symbolically by an apparently random sequence of A^s and B's. Since chaotic trajectories erratically cover 
a large region of phase space, they appear as densely filled regions on a Poincare map. In all three systems 
there is a region of chaos separating the annulus and pretzel regions on the surface of section ^7] . 

A comparison between the relativistic and Newtonian systems reveals differences in the trajectories as 
7] increases. In general, the particles in the relativistic system cross each other at a higher frequency than 
in the Newtonian case for the same value of rj. The structure of the relativistic Poincare maps at all 
energies attainable were similar to the Newtonian ones except for a shifting of all trajectories to one side. 
Remarkably, this structure remained stable up to the values of 77 that were attainable despite the high degree 
of non-linearity in the equations of motion. 

Furthermore, for all trajectories studied in all 3 systems, B motion always occurred in multiples of 3 
jl7| . That is, all symbol sequences were of the form Yli j ^{A™^B^"^ so that any time a single particle, say 
particle 1, crossed the other two in succession, particles 2 and 3 always crossed next before meeting particle 
1 again. 

5.2 Unequal Mass Trajectories 

In order to study the effects of changing the relative masses of the particles, we adopt the parameter a 
whenever two masses are equal so that mi — TO2 — am^. When all three masses are unequal, we will 
describe the relative masses as a ratio (i.e. mi : m2 : 7713 = 1:2:3). 

In the case when two masses are equal {a 1), we find the same diversity of trajectories in the (p. A) 
plane as in the a = 1 case |17j . Figure El shows this diversity for different values of a in the N system while 
Figures and |H1 show examples of annulus and pretzel trajectories respectively for the R system. Similar 
trajectories are obtained when all three masses are unequal. 

In terms of the hex-particle moving in the p-X plane, we did not find a significant difference between the 
equal and unequal mass trajectories besides a general distortion of the annulus orbits as the difference in 
masses increases. For instance, in the equal mass case 117', all of the annulus orbits were generally hexagonal 
about the origin. When the mass of one particle is larger than the rest, these annuli take on a more box-like 
shape, as can be seen in Figured Besides a general distortion, we did not find any novel types of motion 
that were not seen in the equal mass case. Since the qualitative aspects of the motion in the p-X plane does 
not reveal much about the underlying physics, we will forgo any further discussion on this matter. 

Another effect of changing the difference between the particle masses is that the ratio between the number 
of annulus trajectories compared to the number of pretzel orbits at a given energy decreases. That is, as 
a decreases, we find fewer and fewer initial conditions that give annuli compared to initial conditions that 
produce pretzels. The reason for this is best demonstrated by looking at the motion of the particles as a 
function of time. 

Figure ini plots the relative motion of the particles for decreasing values of a in the R system for a specific 
set of initial conditions. We see that, at equal mass {a = 1), a single particle alternately crosses the other 
two without ever crossing the same particle twice, indicative of annulus motion. However, as a decreases, 
the mass of particle 3 decreases and so its frequency of oscillation decreases while its amplitude increases 
with respect to the other two. 

In effect, we see that the two massive particles gravitationally bind together more tightly as the difference 
between their mass and the mass of the third particle increases. Eventually, this binding becomes so tight 
that the two massive particles are forced to execute an additional A motion before crossing the third particle 
and, hence, there is a transition from annulus type motion to pretzel type motion. This behavior, while 
expected for the Newtonian system (Figure Ell j is also present in the relativistic case. 

As the mass difference increases, it is much more difficult to set up initial conditions at a given energy 
such that particle 1 and 2 do not cross more than once during one of particle 3's long period oscillations. 
This effect is also seen in the Newtonian system, as shown in Figure ^| This difference in the ratio of the 
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Figure 6: Examples of trajectories in the N system for different values of rj and a. Each trajectory was run 
for 150 time steps. The small box indicates the starting position of the trajectory. Proceeding clockwise 
from the top left plot the symbol sequences are B, B, AB-^A^B-^, and B^AB^A. 
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Figure 7: Examples of relativistic annulus trajectories for different values of rj and a. Notice the characteristic 
boxy shape at higher values of a. Each trajectory was run for 200 time steps. All trajectories have the symbol 
sequence B. 
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Figure 8: Examples of relativistic pretzel trajectories for different values of rj and a. Each trajectory was 
run for 200 time steps. The top right plot has symbol sequence A^B^. 

number of annulus trajectories compared to pretzel trajectories will be made more clear in the next section 
when we look at the Poincare maps. 

Using these position-time plots, it is interesting to explore the limit where one mass is much greater and 
much smaller than the other two (a 3> 1 and a ^ 1 respectively). The former case is shown in Figure ITTI 
for the N and R systems where a — 100. In both cases, as one would expect, we see the large mass barely 
moves while the other two particles oscillate about it. The inset shows the small perturbations to the motion 
of the larger mass caused by the passing of the two particles. In the Newtonian case, the perturbation is 
very smooth and regular while the perturbation in the relativistic case is more jerky and erratic. That is, 
the velocity of the large mass in the R system increases much more suddenly than in the N system when it 
encounters another, smaller mass. 

Figure El shows the corresponding plots when a — 0.01. We see that, in both cases, the two heavy 
particles form a stable, 2-body sub-system while the third particle oscillates about their center of mass. 
As seen in the upper most insets in both the R and N case, the presence of the light particle has a weak 
gravitational effect, causing the oscillatory motion of the center of mass of the two more massive particles. 
Unlike the effect seen in Figure the perturbation of the motion of the heavy particles due to the crossing 
of the light particle is very small in the R system and almost imperceptible in the N system. The reason 
for this is that the two heavy particles in the a — 0.01 case are twice as massive as the single particle in the 
a = 100 and so the motion of the 2-body subsystem is much more stable and less susceptible to harassment 
from the weaker mass. However, the qualitative nature of the perturbation remains the same as in the 
a = 100 case. 

We also find that the amplitude of oscillations in the Newtonian system is generally larger than in the 
relativistic system at corresponding values of the total energy and that the frequency of oscillations is greater 
in the relativistic case. These observations agree with the results found in the equal mass case (17j . 

Finally, we note that, as in the equal mass case, we find that B motion always comes in multiples of 
three. That is, the symbol sequence always takes the form 

(44) 
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Figure 9: The relative position of each particle with respect to the center of mass is plotted as a function 
of time for various values of a in the R system. The particles 1,2, and 3 have relative masses in the ratio 
1 : 1 : a. Solid line - particle 1, dotted line - particle 2, and dashed line - particle 3. Each plot uses the same 
initial values of (p, X,pp,pX) but the total energy 77 + 1 is fixed by the energy constraint l|12|) . The top two 

plots display annulus motion {B) while the bottom two are classified as pretzel trajectories {{B^AY B^ and 

^253 i^AB'^f respectively). 
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Figure 10: The relative positions of each particle with respect to the center of mass as a function of time in the 
N system. These plots were created using the same procedure as in Figure|^and follow the same conventions 
except that equation H2U|I is used to fix the value of the total energy rj. The first is an annulus trajectory 

(B) while the remaining are pretzels {B^A, {B^A^^B^A^, and {B^A^^f B^A-^ from top to bottom). 
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Figure 11: Relative motion of the particles with respect to the center of mass plotted as a function of time 
for the R system (top) and the N system (bottom). Both plots have mass ratios 1:1:100, or a = 100. The 
lines are as defined in Figure 01 The insets show the small perturbation in the motion of the large mass due 
to the crossing of the smaller masses. 
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Figure 12: Relative motion of the particles for the case where the mass ratio is 1:1:0.01, or a = 0.01 for both 
the R (top) and N (bottom) systems. The insets show the motion of the stable, two body sub-system made 
up of the two heavy particles, as well as the effect of encounters with the light particle. 



for all values of a and 77 that were studied. This extends our hypothesis proposed in Tf, that all trajectories 
in the R and N systems, when translated into a symbol sequence, have the form (|44|l . to also hold for all 
mass ratios of the 3 particles. 

5.3 Global Structure of Phase Space 

By studying the 2 dimensional representations of phase space represented in the Poincare maps we were 
able to discover some interesting global properties of both the N and R systems. We begin this section by 
describing some of the basic features of the Poincare plots and then go on to discuss how the structure of 
phase space changes when the mass ratio of the particles is changed. Our results will then be compared with 
similar studies conducted previously. 

An example of a Poincare map for the N system when all masses are equal is shown in Fieure lTSl All points 
on this surface of section fall within a parabolic region which is defined by the system's energy constraint. It 
was found in |17| . as mentioned previously, that the three types of motion, annulus, pretzel, and chaotic fall 
into 3 regions on the surface of section. Quasi-periodic annulus orbits form single closed loops about a stable 
fixed point at which the motion is completely periodic. In Figure [T^ this region of quasi-periodic annuli is 
located at the center of the plot enclosed within the densely filled triangular shaped region. This densely 
filled region is created by chaotic trajectories and separates the annulus region from the pretzel region. All 
pretzel trajectories fall outside of this chaotic region and form cither a series of disconnected loops or a series 
of disconnected lines. 

A similar segregation of the surface of section is also seen in the R system, an example of which is shown 
in Figure 1141 The relativistic Poincare map is strikingly similar to the Newtonian one in Figure 1131 The 
annulus region is shown as the series of closed loops in the lower right portion, surrounded by the warped 
chaotic region, which is further surrounded by the region of pretzel trajectories. 

In general, the relativistic phase space is a warping of the corresponding Newtonian space. As described 
in jl7| . this is due to the weaker symmetry of the relativistic Hamiltonian compared to the Newtonian. As 



20 





Figure 14: A Poincare map of the relativistic system when the masses of all of the particles are equal. 
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Figure 15: Poincare plots with a = 0.1 for the Newtonian (left) and relativistic (right) systems. The insets 
on the right show the onset of chaos in the pretzel region. 

seen in eq. H2()(l , the N system is invariant under the symmetry pi —> —pi and this is manifest in the symmetry 
about the pn — axis in Figure IT^ The relativistic Hamiltonian, determined by (|12|l . is invariant under the 

symmetry (p,, e) {—pi, — e). Contrary to the Newtonian case, this relativistic symmetry is not manifest in 
our surface of section. 

We find that the annulus and pretzel trajectories continue to fall into similar regions, as described above, 
for all different mass ratios studied, and that these two regions are always separated by a region of chaos. 
By changing the mass ratio at a given value of the total energy, the size and shape of the different regions 
change. 

More specifically, by looking at the case where particles 1 and 2 share the same mass, the annulus region 
becomes smaller and moves towards the top of the allowed region of the surface of section with decreasing 
a < 1 as can be seen in Figure^jfor a — 0.1 in both the N and R systems. This shrinking of the annulus 
region is a manifestation of the effect discussed in t l5.2l whcre annulus motion becomes more difficult to attain 
when one particle is significantly less massive than the other two. 

As a increases, the annulus region extends towards the lower region of the plot, as shown in Figure lTCI for 
a = 10 in both the R and N systems. Essentially, this means that one requires a lower magnitude of angular 
momentum of the hex-particle to attain an annulus orbit in the p-X plane. Since the gravitational attraction 
between the two light particles is not very strong compared to their interaction with the heavy particle, the 
light particles do not tend to oscillate about each other very much but instead act like 2 separate 2 body 
systems with the heavy particle taking the role of the second body, like a two-planet, one-dimensional solar 
system. This situation is shown in Figure ITTl for a = 100 and explains why there is no decrease in size of 
the annulus region with increasing a. 

The symmetry about pr ^ present in Figures [TSlandlTCl is really just an artifact of our choice of surface 
of section. Recall that we chose to construct our Poincare maps by plotting a point each time the hex-particle 
crossed the p — bisector, or, equivalently, each time particles 1 and 2 crossed. The above figures were 
constructed with mi — m2 and so the symmetries of the equal mass system persist. If we were to have 
chosen a different bisector, all of the features discussed above would remain {e.g. shrinking, expanding of 
annulus region) but these would not occur in the same sense and the plots would not be as symmetrical. 

This can be shown by creating Poincare maps for the case when all three masses are unequal. An example 
for both the R and N system is shown in Figure [T7I where the mass ratio is •mi:m2'.ra^ = 1:5:10. Here we 
see that the symmetry about the pn — axis no longer exists in the Newtonian system due to the fact that 
none of the particles have equal mass. We also see a further warping of the relativistic plots due to this 
added asymmetry. Furthermore, we find that the different regions are not as clearly segregated as in the 
mi = TO2 phase space but extend over more of the Poincare map. For instance, in the relativistic map of 
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Figure 16: Poincare plots with a = 10 for the Newtonian (left) and relativistic (right) systems. The insets 
show additional regions of chaos in the pretzel region that are not present in the corresponding region on 
the equal mass Poincare section. 




Figure 17: Poincare plots with a mass ratio of 1:5:10 for the Newtonian (left) and relativistic (right) systems. 
On the left, a) marks the region of chaos separating annulus trajectories (inside) and predominantly pretzel 
trajectories (outside) while the densely filled area directly above and below b) marks a new region of chaos 
amongst the pretzel trajectories. On the right, the densely filled regions marked by a 1) were created by a 
single trajectory separating the annulus and pretzel orbits while the chaotic regions marked by a 2) were 
created by a trajectory within the pretzel region. 
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Figure IT7I we see that the chaotic region separating the annulus and pretzel trajectories (marked by a 1) is 
no longer a single, densely filled loop but actually two loops which were created by a single trajectory. The 
annulus region is confined to the area inside both of these loops, where a single annulus trajectory will visit 
both regions. 

Besides this novel partitioning of the different regions, the changes to the structure of the phase space for 
different ratios of the mass when all three masses are unequal are analogous to the results described above 
for the case where mi = 7712. For example, the ratio 1:5:10 exhibits similar behaviour as the case when there 
is one light particle and two heavy ones, only the effects are not as prevalent because of the intermediate 
mass particle. The results can be seen as an interpolation between the 1:1:1 case and the 1:10:10 case. 

One major difference that we find between the equal and unequal mass cases is the presence of additional 
chaotic regions in the unequal mass space that are not present in the corresponding constant energy hyper- 
surface of the equal mass space. This is true for both the Newtonian and relativistic systems. For the mass 
ratios and energy levels that we have studied, these additional chaotic regions appear within the pretzel 
regions of the corresponding equal mass surface of section. The novel chaotic trajectories are characterized 
by broadened lines in the pretzel region as can be seen in Figures I15H17I The exact, physical mechanism 
that gives rise to this increase in chaos is not presently known. 

Note that, although there are no apparent regions of chaos in the Newtonian, a — 0.1 Poincare map of 
Figure El (besides the one separating the annuli from the pretzels), we do find a slight broadening of lines 
in the other two Poincare maps constructed by plotting points each time the hex particle crosses the other 
two bisectors respectively. 

That is, additional regions of chaos do form in the unequal mass phase space but these new regions 
cannot be seen on the particular choice of Poincare section shown in Figure El We suspect that these new 
regions of chaos would become more prevalent as the difference in particle masses increases. 

It is instructive to compare our results with a similar study of a billiard in i?^ colliding with a wedge 
in a uniform (Newtonian) gravitational field performed by Lehtihet and Miller (referred to as LM herein) 
LM showed that the two-dimensional wedge billiard system is isomorphic to a system of 3 elastically 
colliding, self-gravitating particles (under Newtonian gravity) in one-dimension, with the relative masses of 
the particles directly related to the wedge angle by 



tang=^^ + ^";' (45) 

where a is as defined in our study. LM only considered the situations where the wedge is symmetric, which 
corresponds to the case when two of the three masses are the same. The value oi Q = 7r/6 corresponds 
to a = 1, the equal mass case. This connection between particle masses and the wedge angle agrees with 
the distortion of the potential energy described in § O where the angle of the wedge is related to the angle 
between the bisectors of the hexagonal well. 

The only difference between the LM system and our N system is the existence of collisions in the former 
while the particles pass through each other in the latter. For the case where all three particles are identical, it 
is irrelevent whether one considers that the particles are colliding or passing through each other (besides the 
question of labelling the particles). For this reason, the phase space structure of the Newtonian equal mass 
configuration, as presented in the Poincare maps, is identical between our system and the wedge-billiard 
system. 

LM found that this wedge-billiard system exhibits the characteristics of a conservative Hamiltonian 
system with two degrees of freedom and a discontinuity. By changing the value of a single continuous 
parameter, 0, they found a variety of dynamics similar to our study. More specifically, for 9 < n/A (which 
corresponds to the entire range of physical values of a) they found that integrable, near-integrable (KAM), 
and chaotic regions coexisted in phase space. Furthermore, as the wedge angle was increased from 7r/6 
(corresponding to both an increase or a decrease of a due to the nature of the connection between mass 
ratio and wedge angle (|45|l ). they found that the region surrounding periodic fixed points was consumed by 
regions of simply connected chaos which increased in size with increasing wedge angle. 

As noted above, we find a similar behaviour in both our Newtonian and relativistic systems in that we 
see an increase in the amount of chaos as the difference in the masses increases. However, we have only 
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studied moderate particle mass differences in order to characterize the general nature of the unequal mass 
system and it is not clear how the global structure of our system will behave for very large differences in the 
particle masses. In particular, we do not know if our systems will experience a global transition to chaos 
or if there exists integrable and near integrable regions for all mass ratios. This remains an area for further 
study. 

6 Discussion 

We have presented the results of a continued study of the 3-body problem in lineal gravity begun in 1 161 117| . 
The focus of the present investigation was to see what happens to the motion of the particles when the 
relative masses of each are not equal. Here we summarize our results. 

The derivation of the 3-body Hamiltonian by canonical reduction of the action was summarized and 
the associated post-Newtonian and Newtonian Hamiltonians presented. Each Hamiltonian possesses two 
spatial degrees of freedom with two conjugate momentum degrees of freedom and these were made manifest 
by changing to (p, A) coordinates. Expressions for the potential energy of each system were derived and the 
distortion of the potential energy due to varying the mass ratio was described. 

The results of the study of the equal mass case were summarized and the different types of motion were 
classified into three categories: annulus, where each particle always crosses the other two in succession; 
pretzel, in which two particles can cross each other twice in a row; and chaotic, where the sequence of 
particles crossings does not progress in a discernible pattern. By studying the motion of the 3 particles and 
their corresponding hex-particle representation in the p-X plane, we characterized how changing the mass 
ratio of the particles effects the dynamics of the system. More specifically, we described in physical terms 
how the type of motion (annulus, pretzel, and chaotic) and their relative abundance in phase space changes 
with respect to the mass ratio. 

As the relative difference between the masses of the particles increases, we find the onset of additional 
regions in phase space of chaos that are not present in the equal mass system — in other words, we find 
that motion that was once quasi-periodic is now chaotic. This shows that the unequal mass phase space is 
not simply a deformation of the corresponding equal mass space but, indeed, contains novel dynamics. The 
physical mechanism behind this phenomenon is currently unknown. 

This is similar to the behaviour of a billiard colliding with a wedge (which is isomorphic to 3 particles 
elastically colliding on a line under their mutual, Newtonian attraction) studied by Lehtihet and Miller [2]. 
It is still not known what happens to these novel regions of chaos as the difference in mass gets exceedingly 
large. 

There are still many open areas of study in the lineal, 3-body problem. As was done in the 2-body problem 
|HliniE3^2 it will be interesting to see the effect of adding charge to the particles and a cosmological constant 
to the system. Furthermore, more sophisticated numerical techniques need to be introduced in order to probe 
the dynamics of the system at high energies and in order to study the motion in the post-Newtonian system 
for unequal mass particles. A description of the global structure of phase space for extreme differences in 
the particle masses is still needed in order to determine the stability of the system in these limiting cases. 
As mentioned above, a discrete map between particle crossings in the N and R systems (although it is 
doubtful whether this can be obtained for the latter) may illuminate some of the more general features of 
the 3-body system. The development of a relativistic 3-body system where the particles elastically collide 
instead of passing through each other would also be an interesting subject to study to see if the increased 
chaos reported in 2 has an analogue in the relativistic system. 
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